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Abstract 

In this article we study the effect of a (^'-interaction on a polymerized membrane of arbi- 
trary internal dimension D. Depending on the dimensionality of membrane and embedding 
space, different physical scenarios are observed. We emphasize on the difference of poly- 
mers from membranes. For the latter, non-trivial contributions appear at the 2-loop level. 
We also exploit a "massive scheme" inspired by calculations in fixed dimensions for scalar 
field theories. Despite the fact that these calculations are only amenable numerically, we 
found that in the limit of D ^ 2 each diagram can be evaluated analytically. This property 
extends in fact to any order in perturbation theory, allowing for a summation of all orders. 
This is a novel and quite surprising result. Finally, an attempt to go beyond D = 2 is 
presented. Applications to the case of self-avoiding membranes are mentioned. 
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1 Introduction 



Interacting lines and more generally manifolds play an important role in many areas of modem 
physics. Examples of lines are self-avoiding polymers [1], vortex-lines in super-conductors [2], 
directed polymers in a disordered environment, also equivalent to surface growth [3], diffusion 
of particles and many more. Generalizing results to membranes often poses severe problems, 
but also new insight into physics. Recently, a lot of work has been devoted to self-avoiding 
membranes (see [4] for a review). Applications reach as far as high-energy physics, where 
strings and M-branes have been proposed as a general framework for unifying all fundamental 
interactions. 

In this article, we focus on the statistical physics of systems without disorder. We think 
of situations like the binding and unbinding of a long chain as e.g. a polymer or a membrane 
on a wall or the wetting of an interface. More precisely, we study the interaction of a single 
freely fluctuating manifold with another non-fluctuating, fixed object. Depending on whether 
the interaction is attractive or repulsive, one can distinguish two different scenarios: One may 
either observe a delocalization transition from an attractive substrate as in wetting phenomena 
or steric repulsions by a fluctuating manifold. Both cases have in common that excluded volume 
effects become important. The situation for polymers or a 1 -dimensional interface is relatively 
simple: One knows that in this case polymers interacting with a defect or a wall (excluding 
or not half of the space) as well as short-range critical wetting are in the same universality 
class [5,6]. For two-dimensional interfaces, the situation is more complicated, and a lot of 
effort has been spent to understand e.g. the delocalization transition [7-14]. Particularly, one is 
interested in critical wetting for the case of a short-ranged interaction potential. Then, it can not 
be approximated by a polynomial in the field r and the conventional field-theoretic approach 
known from 0"^— theory fails. This led to a couple of different ansatze [8,9,11-13], among 
others the functional renormalization group approach of [12,13]. 

Here we follow a different route: We start by constructing a systematic perturbative renor- 
malization group treatment of the delocalization transition as well as the universal repulsive 
force exerted by a membrane on a point, line or more generally hyper-plane like defect. We 
therefore introduce a flexible "phantom" manifold of internal dimension < D < 2. By intro- 
ducing a 5-potential in a subspace S of dimension d' , part of the embedding space W^^'^ (see 
fig. 1.1) we punish configurations crossing 8. Neglecting the effect of self- avoidance between 




Figure 1 . 1 : A fluctuating membrane interacting by excluded volume with a surface-like defect. 
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distinct points within the fluctuating manifold, the free energy of a given configuration reads: 

n[T\ = [ d^x kvr{x)f + go [ d^x [ d^'y 6'^'' {r{x)-g) , (1.1) 

where f{x) is the position of monomer x G in embedding space W^'^'^' , and Ai denotes 
the D-dimensional coordinate space of the manifold. In the case of a polymer this is simply 
the internal chain length, qq is the attractive or repulsive interaction, of dimension (in inverse 
length-units) 

e{D,d):=[go] = D-'^-^d. (1.2) 

The interaction is naively relevant in the infrared for e positive, and irrelevant for e negative. 
In a perturbative treatment, e is the natural expansion parameter. The situation is similar to 
self-avoiding membranes and thus (1.1) has been studied as a "toy-model" for the analysis of 
the renormalizability of the more complicated interaction in the generalized Edwards model for 
self-avoiding polymerized membranes [15-17]. It was shown to be renormalizable for arbitrary 
manifold dimensions 0<D<2 [16,17]. This means that the large scale properties are universal, 
i.e. do neither depend on the regularization scheme used in these calculations, nor on the form 
of the contact interaction, as long as it is short-ranged. Universal quantities have thus been ob- 
tained to one-loop order. They are related, as we will show below, to the correction to scaling 
exponent uj, which as usually in critical phenomena contains deviations from the long-distance 
scaling behavior. Neglecting self-interactions, one has to distinguish between polymers, which 
are always crumpled and membranes, for which a highly folded high-temperature state is sep- 
arated through a crumpling transition from an essentially flat low-temperature phase [18-21]. 
The effect of self-avoidance can also be taken into account [22-25,4]. It has been shown to lead 
as in the case of polymers to a swelling of the membrane and correspondingly to a non-trivial 
exponent for the radius of gyration [26,27]. However, it is not clear whether a fractal phase can 
be found in experiments or simulations. (See [28] for the latest simulations). In this work we 
will mostly neglect self-avoidance. 

The aim of this article is two-fold: First, we present the necessary techniques to treat the 
model (1.1) beyond the leading order. We explicitly perform a two-loop calculation which gives 
the correction to scaling exponent uj at order e^. Specializing to membranes one finds that the 
2-loop result naively diverges in the limit oi D 2. This is a problem of the ^-expansion, 
since there, diagrams have to be evaluated at e{D, dc{D)) = 0, and taking D — > 2 implies 
dc{D) oo causing the result to diverge. This motivated us to try a "massive scheme" in fixed 
dimension, i.e. finite e. It turns out that the limit D ^ 2 can then be taken and is regular. Even 
more, the 2-loop diagram, which for D < 2 can only be calculated numerically, can now be 
evaluated analytically. This striking property even holds for higher orders, and we are able to 
give an explicit - quite simple - expression for the perturbation series. Then the whole series 
can be summed and the strong coupling limit analyzed. This is one of the very few cases where 
one can indeed obtain the exact relation between bare and renormalized coupling in the limit of 
D ^ 2, and thus the exact /i-function. This result does not depend on the explicit regularization 
and renormalization prescriptions, and is also obtained for a membrane of spherical or toroidal 
topology. In a final step, we lay the foundations for an expansion about D = 2. In contrast 
to the leading order, the first order corrections already depend on the cut-off procedure. We 
study one specific procedure, which turns out to reproduce results for polymers at leading order 
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approximately, and even exactly in d — 0. Work is in progress to obtain a more systematic 
expansion about D — 2 [29]. 



2 Model and physical observables 



2.1 The model 

We start from the manifold Hamiltonian (1.1). We split the total embedding space M'^+'^' 
into £ and its orthogonal complement £± of dimension d. Each x E V points to a point 
f{x)={f±{x),f\\{x)), with f\\{x) e £ and f±{x) G £±. The integration over the subspace £ 
is then trivial and gives 



6''+'''{r{x)-y) = 6%f±{x)) . 



In the partition function 



^totai = Jvlrl exp(-7i[f]) 
the contributions from the parallel and orthogonal components of f{x) factorize as 

with 



^0 = J V[f\^] exp(^-^{Vf{x)y 



V[f^] exp(-7i^[rl]) 



^±[rj = dPx Q(Vr-l(x))2 + ^0 nr±{x))^ ■ 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



Since Zq is trivial, we will only consider H_l [r±\ and shall drop the subscript j_ for notational 
simplicity. We keep in mind that cases with d < D make sense, for instance a polymerized 
(non- self avoiding) membrane interacting with a wall is described by (2.6) setting D = 2 and 
d^l. 

Let us discuss (2.6) in more detail: The first term is the elastic energy of the manifold 
which is entropic in origin. We have scaled elasticity and temperature to unity. The second 
term models the interaction of the manifold with a single point at the origin in the external 
(i-dimensional space, but we remind that the physical interpretation may well be that of a line 
or surface. The coupling constant go may either be positive (repulsive interaction) or negative 
(attractive interaction). We now give the dimensional analysis. In coordinate-space units, the 
engineering dimensions are 



dim[a;] 
v := dim[r] 



1 

2-D 



e :-- 



dim 



d^x 6\f{x)) 



M 



D-vd 



(2.7) 
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D 




Figure 2.1: Critical line defined through e = <^ dc{D) = ^Bd- ^he interaction is relevant for points that lie 
above that line. 



dim[//] — —1 , 



where 



(2.8) 



is an inverse length scale. The interaction is naively relevant for £ > 0, i.e. d < dc with (see 
figure 2.1) 

dc = , (2.9) 

irrelevant for £ < and marginal for £ = 0. It has been shown [16,17] that the model is 
renormalizable for < D < 2 and e > 0. Results for negative e are obtained via analytical 
continuation. One can define the renormalized coupling g as 



9 — 



V 



[Z(0) - Z{g^)] V , 



(2.10) 



M 



where the normalization J\f depends on the definition of the path-integral (but not on L) and is 
chosen such that 

g^goL' + Oig^,). (2.11) 
Universal quantities are obtained at fixed-points of the /3-f unction, which is defined as 

dg 



Pig) 



(2.12) 



90 



The /3-function thus describes, how the effective coupling g changes under scale transforma- 
tions, while keeping the bare coupling go fixed. Let us already anticipate the 1-loop result. 



H. Pinnow and K.J. Wiese, Interacting Crumpled Manifolds 



7 



which we derive later. It reads 



P{g) = -eg + 1/ + 0{g') , 



(2.13) 



where g is the dimensionless renormalized coupling. Apart from the trivial solution, g = 0, 
the flow equation given by (2.12) and (2.13) has a non-trivial fixed point at the zero of the 
/^-function 



We shall show below that the scaling behavior is described by the slope of the RG-function 
at the fixed point, which is universal as a consequence of renormalizabilty. The long-distance 
behavior is then governed by the (^-interaction as considered in our model (2.6), which is the 
most relevant operator at large scales. Let us now discuss possible physical situations (see fig. 
2.2): 

(a) £>0: The RG-flow has an infrared stable fixed point at g'* > and an IR-unstable fixed 
point at g — 0. The latter describes an unbinding transition whose critical properties 
are given by the non-interacting system, while the non-trivial IR stable fixed point deter- 
mines the long-distance properties of the delocalized state, the long-range repulsive force 
exerted by the fluctuating manifold on the origin - which we remind may be a point, a 
line or a plane. 

(b) £<0: Now, the long-distance behavior is Gaussian, while the unbinding transition occurs 
at some finite value of the attractive potential, 5'*<0, which corresponds to an infrared 
unstable fixed point of the /3-function. Below g* the manifold stays always attracted. 

(c) £ = 0: This is the marginal situation, where the transition takes place at g* — 0; we 
expect logarithmic corrections to scaling. 

Note that in the presence of an impenetrable wall constraining the configurational space strictly 
to half of the embedding space, the above considerations should still apply, when shifting the 
interaction strength appropriately. We shall discuss that in 2.3. 



g* = 28 + 0(6'). 



(2.14) 





g 



g 







g 



(a) 



(b) 



(c) 



Figure 2.2: RG-function and flow for increasing manifold size L for the dimensionless renormalized coupling g: 
(a) in the case £>0, (b) in the case s<0, (c) in the case £ = 0. 
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2.2 Repulsive force exerted by a membrane on a wall 



Since we are mainly interested in the long-distance properties of membranes for which always 
£ > (this is (a) on fig. 2.2), let us try to calculate the repulsive force exerted by the membrane 
on the origin in the case that this point is strictly forbidden. We will derive a universal expression 
for this force [17]. We need the (not normalized) membrane density at position r 

Z{f, 9o) - / ^[^ / S^'ifix) - r) exp[-H] . (2.15) 

J Jm 

Since the (5-interaction also appears in H, we can relate the density at the origin to the derivative 
of the partition function with respect to qq: 



1 d 



Z[go) = -r? — 
M dgo 



(2.16) 



where g is the renormalized or effective coupling defined in (2.10). Since Z{0, go) is dimen- 
sionless, it depends on the dimensionless combination goL^ (and e) only. Using /j, — 1/L, we 
obtain 



AfZ{0,go) = L-^^ 
dgo 



(2.17) 



For the further considerations, it is sufficient to know the behavior of the /3-function close to the 
nontrivial zero. Expanding about g—g*, we have 



m = (9 - 9*Mg*) + 0{{g - g*f) . 
Combining (2.16), (2.17) and (2.18), we get 



dg 



egoL^ 



{g - g*Mg*) 

sgoL^ 



d(goL-) 

The solution of this differential equation is 

g-g*^ {goL^)-^^^^^'' 
Using (2.17), this implies the scaling law 

Z{Q,go)r^{gl''L)-^^'*^-^ . 



+ 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



We now use this result to derive a scaling law for Z{r, go), with f ^ 0. In addition to goL^ 
(and s), Z(f, go) depends on r; taking into account rotational symmetry and dimensionality, it 
depends on r only through the dimensionless combination r/L'', with r = |r|: 



Z{f,go)^Z{r/L^,goL' 



(2.22) 



For Z{f,go), the most interesting limit is that of g'o — • Physically, this corresponds to strictly 
forbidding monomers to be at the origin. Therefore it is clear that this limit is well-behaved, 
and 

Zoo (r/L^) := lim Z(r/L^yo^^) (2.23) 

go-*oo 
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is finite. In order to be consistent with (2.21) it has to obey a power law in the scaling regime 

Z^{T/r) ~ {r/L^f . {12A) 
Comparing the L-dependence of Z^{r/L^) and Z(0, qq), we obtain the exponent identity 

e = . (2.25) 

Finally, from (2.24) we derive the repulsive force between the origin and the manifold 

— * 

/(r) = Vr\nZ^{\r\/U) = 9-^ . {1.16) 

Note that to derive this result, ksT has been set to 1. Reestablishing the temperature-dependence, 
we find 

— * 

/(f) = ksTe^ . {Ill) 

Also note that this argument gives ^ = at the Gaussian fixed point, which is necessary since 
for g^o = no force is excerted on the membrane. 

2.3 Unbinding transition 

Let us discuss the physical situation at the UV-stable fixed point in figure 2.2. The fixed point 
corresponds to a delocalization transition of the manifold, which is at vanishing coupling g*—0 
for e>0 and at some finite attractive coupling g*<0 for €<0. 

In the localized phase g<g*, correlation functions such as {[f(x) — f(y)]^) and the associated 
correlation length (in the D-dimensional internal space) should be finite, as well as the radius 
of gyration ^j^. Approaching the transition point these quantities diverge as [14] 

Cii ~ (g^-g)-"^^ , U ~ (g^-g)-"^ ■ (2.28) 

Since ~ ^||', the exponents z/y and iy± are related through 

= , (2.29) 

u being the dimension of the field (2.7). To relate z/y to cu, we first observe that ^\\ depends 
only on go and the membrane size L, or its inverse n = 1/L. Writing 4|| as a function of the 
renormalized coupling g and the scale we thus have 



-{-f{g) + mf{g)) , (2.30) 

A* 
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Figure 2.3: A polymer (here for simplicity directed) in interaction with a wall (thick dots). Left before and right 
after flapping up. The grey area is impenetrable. 

where the dimension-full factor 1 //x has been factored from the dependence on the dimension- 
less renormalized coupling g. Using that ^{|| (g) = j[f'{g), this can be rewritten as 

^ u;{g*){g - g*)^^\^{g) . (2.31) 

The solution to the above equation is 

eil-l^-Zr^"^"*^- (2.32) 

This leads to the identification 

^\\ = — , ^± = — • (2.33) 

^[9 ) ^[9 ) 

Note that a;(g'*)<0 at the transition. Specializing to (D, ci)=(l, 1), we find 

1/^ = 1 , i^ii = 2 . (2.34) 

These exponents are also valid for the delocalization transition of a 1 -dimensional interface 

from an attractive hard wall in 2-dimensional bulk space [8,14,6]. This can be understood 
as follows: The partition function of a fluctuating polymer interacting with a (5-defect at the 
boundary of a hard wall can be written as 



lf:^4lT ^ (2.35) 



2^ V2 

n=0 ^ 

where the sum runs over the number of contacts with the defect hyper-plane. E denotes the 
contact energy, p = {kBT)^^ and Zn is the constrained partition function of the free interface 
having exactly n points of contact with the defect. The powers of ^ only appear in the presence 
of an impenetrable defect and reflect the fact that all configurations in the presence of a hard wall 
can be obtained by flapping up those parts of the polymer which have penetrated the line (see 
figure 2.3). This accounts for a factor | for each possible flap. Thus, the hard wall constraint 
translates into a shift in the binding energy according to 

Eno wall -Swall — E^o wall + ksT In 2 . (2.36) 

The delocalization transition occurring for {D, d)={l, 1) at vanishing potential strength is shifted 
to an attractive interaction dependent on the temperature. Due to the correspondence (2.36) the 
exponents characterizing the transition remain unchanged. 



H. Pinnow and K.J. Wiese, Interacting Crumpled Manifolds 



11 



3 Operator product expansion and one-loop result 

The partition function of the model has been defined in (2.5). Analogously, we define the 
expectation value of an observable 0{r) as (we denote r for r from now on) 

iOirK ■= ^^"-l . (3.1) 

Most physical observables can be derived from expectation values of general M-point vertex 
operators like 

M 

0{r) = Y[e'P^'^'^K (3.2) 

Before we go on, let us make some changes in normalizations which will be helpful in the 
following. First, we rescale the field and the coupling constant according to 

r^{2-DY/^r, ^ . (3.3) 

Second, we change the integration over internal coordinates to 



1 r 7r^/2 



being the volume of the D-dimensional unit-sphere. Third, the normalization of the ^-distribution 
is changed to 

5'^{r{x)) := (47r)'^/25(r(x)) = / e^'^^^^) (3.5) 

Jk 

with 

= Tf-'^/^ / d'^k . (3.6) 



The advantage of these normalizations is that 

e-'^' = 1 , / e(|a;| < 1) = 1 . (3.7) 



Ik 



'k 

The model in the new normalizations now reads 



n[T] = ^ \ {VT{x)f + goj^ ~S'{r{x)) , (3.8) 
and (due to the factor of 1/(2 — D) in the above definition) the two-point correlator is 

C(x,-xj) := ^ {[r{x,) - rix^)]')^ = \x, - x^^ ■ (3-9) 

We now proceed with the calculation of physical observables. As an explicit example, let 
us consider the perturbation expansion of the 1-point vertex operator 

Z^) := Z{g,){e^^<^)),, = J] ^-j^Z^, , (3.10) 

N=0 
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where 

This can be written as 

N 



zip) 
■'N 



i=l "'^i ; 



(3.11) 



i=l 



N 



(^kir{xi) +pr{s)) 



(3.12) 



where we have already integrated out a global translation of the field r. The Gaussian average 
is: 



cxp 



N 



i(^kir{xi) +pr{s)) 



1=1 



cxp 



1 ^ 

9 kikjC{xi—Xj) 



, ko^p. 



(3.13) 



The Z^^ in (3.11) posses short distance singularities for N>2. In order to analyze these sin- 
gularities we will use the techniques of normal ordering and operator product expansion in the 
sequel. In our problem the procedure of normal ordering a vertex operator e^'^^^^^^^ turns out to 
be quite simple: In any expectation value, we factorize out the contractions between the vertex 
operators. At one-loop order or equivalently at second order in (3.10) this means that 



Qikir{xi}^ik2r{x2} _ . Qi{kir(xi)+k2r(x2}) . Qkik2C(xi-X2} 



(3.14) 



This can be understood as a definition of the normal ordered product :e^i'<^i^(^i)+'<^M^2)) . Note, 
that equality in the above equation is meant in the sense of operators, that is when inserted into 
expectation values within the free theory. 

Let us now turn to an explicit derivation of the OPE of two 5-interactions, i.e. we study (3.14) 
for small distances X1—X2. Since the short-distance singularities appear only in the internal con- 
tractions, the normal ordered product in the r.h.s. of (3.14) is regular for small distances X1—X2 
and thus can be expanded therein. We then project the resulting terms on the corresponding 
operators. For that purpose we make a change in the momentum variables according to 



'^ + k2 



(3.15) 



In internal space we change coordinates to the center of mass system: 

X := -{xi + X2) , y:=X2-xi. 
Thus, the r.h.s. of (3.14) takes the form 



(3.16) 



,[(^_fc2)r(S-3//2)+(^+fc2)r{S+2,/2)] . (kj /4-k^^)C{y) 



.^ikir{x) 



1 ik2yS/r{x) i/ci — — — VaV/3r(x) 



{k2yVr{x)f 

2r 



1 + iC{y) 



2! 4 



+ 



X 



(3.17) 
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Integration over ki, k2 yields the OPE of 5'^{r{xi)) with 5'^{r{x2)): 

~5M^ - |))5~'(r(x + |)) = y-'~5\r{x)) - y-'^'^\ dl^H^)) 

_^-.o!-2.+2_L^d(^(^))(Vr(x))2 

dj\r{x))^n{x) + (3.18) 

Restricting ourselves to the most relevant next to leading operators for internal dimension D 
smaller than, but close to 2, the operator product expansion reads: 

¥{r{x - l))¥{r{x + |)) = y-^'SM^)) + ^-'-^+^(^-2) (-g.^mr(^)) _ ^3 ^9^ 

n=l 

The OPE-coefficient in front of the leading operator, which is the 5-interaction itself, carries the 
leading short distance singularity. In the following we introduce some diagrammatic notation. 
Denoting • := 5'^{r{x)) we abbreviate the projection of two nearby (5-interactions on the 5- 
interaction itself as 

((• •)!•) = y-^". (3.20) 

We now analyze short-distance divergences of the perturbation expansion using the OPE in 
(3.18). Being interested in the leading scaling behavior, we drop all subdominant terms for 
< D < 2, such that we only need the leading OPE-coefficient (3.20). The integral over the 
relative distance in ( • • |») • is logarithmically divergent for £=0. In order to define a counter 
term we use dimensional regularization, i.e. we set £>0. The integral is then UV-convergent, but 
IR-divergent. An IR regulator or equivalently a subtraction scale iJt,—L~^ has to be introduced to 
define the subtraction operation. Generally, we integrate over all distances bound by the inverse 
subtraction scale L. The OPE-coefficient in (3.20) then reduces to a number as 

(( • •:»^ := / ((• •, !•) = L'f{e, D) . (3.21) 
J\y\<L 

Let us use the scheme of minimal subtraction (MS). The internal dimension of the membrane 
is fixed and (3.21) is expanded in a Laurent-series in e, starting here at order e~^. Denoting the 
term of order £^ in ( | by ( | the simple pole in (3.21) is 

(€«)hl-i = J- (3-22) 

Of course, in this case, due to our normalizations, this is exact. 

Let us rewrite the Hamiltonian in (2.6) in terms of the renormalized dimensionless coupling g. 
For perturbative calculations it is custom to write 



(3.23) 
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where Zg is the renormalization factor, which is fixed by demanding that observables remain 
finite in the limit of £ ^ 0. The Hamiltonian becomes 

We find to one-loop order 

= 1 + y^. .^^., + 0{g^) . (3.25) 

There is no field-renormalization since the elastic part has also to describe the manifold far 
away from the origin, for which the interaction term can be neglected. Formally, of course, this 
follows from the OPE. The renormalization group /3-function is defined as 



9 , (3.26) 

so 



from what follows together with (3.23) that 



-£9 



13(g) = 1^ . (3.27) 

To one-loop order we obtain from (3.27) and (3.25) as anticipated earlier 

Pig) = -eg + ^ + 0{g') . (3.28) 

The long-distance behavior of the theory is governed by the IR-stable fixed point of the RG- 
flow. Fixed points are zeros of the /3-function: 

(3{g*) = g*=0 or g* ^ 2e . (3.29) 

In section 2.2, we had defined the correction-to-scaling exponent lo. It is obtained from the 
slope-function u;{g), defined as 

u^ig) ■■= ■ (3.30) 

The correction to scaling exponent is u;{g) evaluated at the fixed points, and is found to be 

u;{g*^0)^-s, uj{g*^2e)^e. (3.31) 

Of course, this result is apart from the very existence of a fixed point rather trivial, since it is 
determined through the leading term of the /3 -function, which is always —eg. In order to obtain 
any non-trivial information from the perturbation series it is thus necessary to go beyond the 
leading order. 
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4 Two-loop calculation in a MS-scheme 

4.1 Operator product expansion 

Let us now continue with the calculation at the two-loop order. At each order of perturbation 
theory there is only one new diagram. At two-loop order this comes from three coalescing S- 
interactions. Again, let us rewrite the product of the three vertex operators as its expectation 
value times the normal ordered product: 



^ikir{xi) ^ik2r{x2) Qik3r{x3) 

— . Qi[kir{xi)+k2r{x2)+k3r{x3)] . Qkik2C{xi-X2) Qkik3C{xi-X3) Qk2k3C{x2-X3) 



(4.1) 



We use the following change of coordinates 



X := ^{xi + X2 + xs) 
a:— X2 — Xi 
b:^X3- X2 
c:— Xi — Xs 



Xi = X + ^{c — a) 
^ ^ X2 = X + ^{a - b) , 
X3 = X + ^{b - c) 



(4.2) 



As at 1-loop order, :Q^['^M^i)+k2r{x2)+k3r{x3)] . fj.gg divergences upon approaching the points 
Xi, whereas the UV-divergence comes from the factor e'^ik2C{xi-x2)Qkik3C{xi-x3)Qk2k3C{x2-x3)_ 
The leading term is thus obtained upon setting xi = X2 = .x'3 = x in :c*[^i^(^i)+'^2r(x2)+fc3r-(a;3)] 
dropping the gradient terms. Further shifting ki ^ ki — k2 — k^, we obtain 



gifcir(i) ^{ki-k2-k3)k2C{a)+k2k3C{b)+{ki-k2-k3)k3C{c) 



(4.3) 



Integrating over ki, k2 and /C3 yields the final result 
S''{r{xi))S''{r{x2))S''{r{xs)) = (i/,;|«)^'^(r(x)) + 



(4.4) 



where the OPE-coefficients are given by 



C{a)C{c)--{C{a)+C{c)-C{b)r 



-d/2 



-{a^+b^+c") {b^+c'-a") {a^+c'-b") {a'W-c") 



-d/2 



(4.5) 



which contributes to the renormalization factor at two-loop order. The first subleading term is 
(denoting := {-dl)~6\r{x))) 



-\c{a)C{b)C{c) 



{a'+b'+c") {b^+c'^-a") {a'+c'-b") {a'^+b^-c") 



-d/2-1 



(4.6) 
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4.2 Numerical calculation in < L> < 2 



Let us now turn to the explicit calculation of the second order contribution to the coupling con- 
stant renormalization. To two-loop order the renormalization group Zg-factor for the coupling 
constant reads 



6 



i- 



--( 

-1 2\ 



• • • 



)U+0{g'). 



(4.7) 



Note that this Zc,-factor can either be obtained from (2.10) and (3.23) or by expanding the 
perturbative series and checking that divergent contributions proportional to 5{r{x)) cancel. 

The formula is arranged such that in the brackets the second order pole in e cancels. This is, 
because to leading order in 1/e, the two-loop diagram factorizes into two one-loop diagrams. 
The combinatorial factor of 3/2 is composed of a factor of 3 for the number of possible subdi- 
vergences, and a factor of 1/2 for the nested integration: The subdivergence in the distance say 
a only appears in the sector (i.e. part of the domain of integration) in which a is smallest. 

Let us define: 

3 , , ,2 



h{D,L):-- 



(4.8) 



Let us already note that we will later calculate in a "massive" scheme in fixed dimension d. 
In order to use the results derived here, we keep d arbitrary instead of setting d = dc- We 
now want to derive an expression which can be integrated numerically. Since the bounds on 
the integrations in (4.8) are asymmetric (remind that in ([/,) |«)^ all 3 distances a, b, and c 

are bounded by L, whereas ( • • |»)^ is an integral with two bounds only) we treat both terms 
separately. We start with the first one, which can be written as 

I^{D, L)= [ e(a < L)G(6 < L)Q{c < L) 1 •) . (4.9) 

Ja,h,c ^ ^ 

The divergence when integrating over the global scale is eliminated by noting that L^l2{D, L) - 
2e hiD, L). There remains of course the sub-divergences, which will be subtracted later by the 
counter-term. Applying this procedure to (4.9), we obtain 



2elf{D, L)= I 8{a - L)Q{h < L)Q{c < L) + 5{c - L)Q{b < L)Q{a < L) 

J a,b,c 

+5{b - L)e(a < L)e(c < L) (|(/,) •) . (4.10) 

Since the integral is symmetric in a, b, and c, we can replace it by 3 times the first term. 
Introducing proper normalizations, an explicit representation for the measure and factoring out 
the explicit L-dependence, we obtain 



o /-oo /-oo 

2eI^{D, L) = 3L^'^^ / dc,, / dc^c^-^Qib < l)Q{c < 1) ( /, ; •) 

J-oo Jo ^ ^ 

We also have to derive an expression for the counter-term. It can be written as 



a=l,b,c 



(4.11) 



-'2 



{D, L) = -- ! I a-^V^"^ e(a < L)Q{c < L) . 

JaJc 



(4.12) 
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Figure 4.1: Area of integration (dashed) of the first integral in (4.13). The intersection with B is to be mapped onto 
A 



Using the same trick to derive w.r.t. the IR-regulator as above, we obtain 

3 



2eI^{D,L) = ~L'' 



A,B,C 



6{a - l)e(c < a) + 6{c - l)Q{a < c) 



. (4.13) 



Using the technique of conformal mapping, which is presented in appendix B, both integrals 
above can be transformed into one having support in the single sector A defined through 
h,c < a, (where here a is set to 1). Analogously, we define sectors B and C, as the subset 
of domains where h respectively c is the longest distance. Note that the first term in (4.13) only 
has contributions from A and B, while the second can be restricted to A and C after renaming 
a into h and c into a. 

Performing explicitly the mapping of both B and C onto A one obtains 



/ {a-^^c-''^Q{c < a)) 
Jb 



a=l 



= h 



-vd—2e„—vd 



A 



[ {a-'"^b-'"^e{b < a))ia=i = I c-"^-^'b-'"^Q{b < c) . 

Jc J A 

Combining the contributions from all sectors, we arrive at 

2eI^{D,L)^ --L^' [ {b-'"^+c-'"^+b-'"^c-'"^maK{b,c)-'^') . 
2 Ja 

Combining (4.11) and (4.15), we obtain an expression for the complete diagram 

Q PCO poo 

2el2iD, L) = 3L^'^^ / dcy / dc^c^^^Q^^ < ^^q^^ < 

J-oo Jo 



X 



(4.14) 



(4.15) 



-{l+b''+c''){b''+c''-l){l+c''-b''){l+b''-c'') 
- [6-'^'^+c-'^'^+6-'^'^c-'^'^max(6,c)-'"] | . (4.16) 



This expression will also be used later in a massive scheme. Since here we are only interested 
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in the £-expansion, we take the limit of £ —> in the integrand. We have made the (non- 
trivial) check that the integrand remains a finite integrable function, which can be integrated 
numerically: 

T 2s Q poo i>oo 

h{D, L) = 3—^ / dc|| / dcxc?-2e(6 < l)e(c < 1) 

2.£ bo J-oo Jo 



X 



- (l+6"+c") (ft^+c'^-l) (l+c"-6") (l+ft^-c'^ 



-| -dc/2 



(4.17) 



Let us now proceed with the explicit numeric evaluation of (4.17). There are integrable sin- 
gularities both in C|| and c±. To disentangle them [30,27] we make a change of variables from 
Cartesian coordinates (cy, c_l) to radial and angular coordinates (c, a) according to 



CCDS q; 



c_L = csmlo; ) 



(4.18) 



such that 



/oo rc 
dC|| / 
-oo -'0 



dc±c 



D-2 



1-1 /■7r/2 

/ dcc-^-^ / da (sin(a)) 



D-2 



(4.19) 



where the upper bounds of 7r/2 on a and 1 on c is due to the bound of 1 on 6 and c. We can 
further restrict integration to the half-sector with c < 6. In this sector, there remain singularities 
for small a and small c. The reader is invited to verify that they are eliminated by a second 
change of variables 



a — 



1 




D 




D 




D 




0.1 


6.71 


0.8 


-0.102 


1.5 


2.24 


0.2 


2.11 


0.9 


-0.0712 


1.6 


4.45 


0.3 


0.847 


1.0 





1.7 


10.4 


0.4 


0.325 


1.1 


0.123 


1.8 


36.7 


0.5 


0.0807 


1.2 


0.321 


1.9 


658 


0.6 


-0.0415 


1.3 


0.644 


2.0 


oo 


0.7 


-0.0949 


1.4 


1.20 







Figure 4.2: Numerically obtained results for the two-loop diagram in the MS-scheme. Some values are tabulated. 
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(D,d)=(l,l) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
D D D 



Figure 4.3: Correction to scaling exponent to for {D=l,d=l) with different extrapolation parameters 
{D, rf), {D, e{D, d)), {D, Dc{d)) (from left to right). The abscissa labels the corresponding internal dimension 
Do of the departure point on the critical hne e{Do, do)=0. 



The above method works for D > 1 only. For D < 1, one has to analytically continue the inte- 
gration over c± in (4.11). This can be done by partial integration, since the integrand does not 
explicitly depend on c± but on c = ^c]_ + c|. To have no boundary terms, one should partially 

integrate before mapping everything onto the sector A. We have explicitly checked consistency 
with the earlier used measure for 1<D<2. Note finally, that in the special case of D = 1 in 
(4.17) the integration measure becomes a distribution reducing the integration measure to a line 
along C||. 

The diagram is shown on figure (4.2). It grows as for D ^ and exponentially for 

HD)r.i^-j r^l^-j . (4.20) 

In the case of polymers D=l it vanishes. The reason is that the diagrams factorize, and thus the 
1-loop result is correct to all orders when working in an appropriate scheme. (See e.g. [4].) 

4.3 RG-function and extrapolation 

^From (3.27) the renormalization /^-function in second order of perturbation theory reads 

2 3 

Pig) ^-eg + ^- ^h{D) + 0{g') . (4.21) 
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(D,d)=(2,l) (D,d)=(2,2) (D,d)=(2,3) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



(D,d)=(2,l) (D,d)=(2,2) (D,d)=(2,3) 




0.2 0.4 0.6 0.8 ^ 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 ^ 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 ^ 1.2 1.4 1.6 1.8 



{x,y) ; {D,d) 


(1,1) 


(2,1) 


(2,2) 


(2,3) 


{D, Dc) , max 


0.55 


2.80 


2.43 


2.10 


{D,d) , max 


0.55 


2.80 


2.43 


2.20 


{D, e) , max 


0.55 









Figure 4.4: Correction to scaling exponent oj for different extrapolation points {D,d): Extrapolation parameters 
are {D,Dc{d)) (top) and {D,d) (bottom). The abscissa labels the corresponding internal dimension Dq of the 
departure point on the critical line e{Do, do)=0. Furthermore, the two-loop values with minimal sensitivity with 
respect to the departure point are tabulated. 
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The long-distance behavior of the theory is governed by the IR-stable fixed point of the RG- 
flow. Fixed points are given by zeros of the /3-function: 

P{g*)^0 g*^0 V -e + ^-^l2{D) + O{g*^)^0. (4.22) 

The physically interesting nontrivial fixed point is the one that is next to zero. In an ^-expansion 
it reads: 

g*^2e + ^l2{D)e' + 0{e^) . (4.23) 
The correction to scaling exponent u at this fixed point is found in an ^-expansion to be 

cuig*) = e - p2{D)e'' + 0(8^) . (4.24) 

The question now arises of how to calculate u for a given physical situation. The general 
strategy that we follow in the MS-scheme is that we make use of the freedom to choose any 
point on the critical line (2.2) around which to start our expansion [31,30]. To expand towards 
some physically interesting point {D,d) we furthermore dispose of the freedom to follow any 
path starting on the curve £=0. Since our expansion (4.24) is exact in D but only of second 
order in e, we can equally well expand it to second order in D around any point on the critical 
curve. Now we can change our extrapolation path through an invertible transformation (x, y) — 
{x{D, d),y{D, d)). One expresses D and e as functions of x and y and reexpands uj to second 
order in x and y around the point (xq, yo) — {xo{Do, dc{Do)), yo{Do, dc{Do))), where we recall 
that dc{D) is defined such that e{D, dc{D)) = 0. The aim is to find an optimal choice of 
variables {x,y). The guidelines for such a choice have been discussed in [27], where this 
procedure has been successfully applied to extrapolate results for the anomalous dimension of 
self-avoiding membranes. As a general rule, we "trust" most a result that is the least sensitive 
to the starting point of the extrapolation. Therefore, we would like to find a plateau for some 
suitable chosen variables (x,y). This procedure works well for polymers, or more generally 
for membranes with inner dimension close to 1 . However it turns out that for membranes it 
works much less well than in the self-avoiding case: we could not find a plateau, but at best an 
extremum. 

Some extrapolations are shown in the following figures (4.3, 4.4). We start with polymers 
interacting with a 5-wall corresponding to the point {D, d) — (1, 1). (Note that for d < 2 the 
interaction becomes relevant for polymers.) This can be seen on fig. 4.3. The values obtained 
from the plateaus appearing at the two-loop level (cj?^i0.55) are quite close to the exact value 
UJ = 1/2, known from the fact that the result in D = 1 at one-loop is exact. From the latter it 
follows that it is best to use (D, e) as extrapolation parameters and to stay in D — 1. 

Let us now turn to the results obtained for uj extrapolated to points D — 2 (fig. 4.13): The 
extrapolated value depends strongly on the starting-point on the critical line. We can identify 
a least sensitive value, which lies below D = 1, reflecting the D-dependence of the two-loop 
diagram. We obtain identical results for different extrapolation paths. Note, that if we fix d = 1 
and perform an £-expansion in D as is commonly done [32], then, this corresponds to starting 
from Dc{d—1) — 2/3 < 1. The d-dependence of the extrapolated values is recovered in the 
massive scheme, which we discuss in the next section. 
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5 Calculation in fixed dimension: The massive scheme 

5.1 Phantom manifolds 

An alternative scheme to an ^-expansion is a calculation in fixed dimension. This way we hope 
to circumvent the difficulty that we faced when calculating the two-loop diagram in the previous 

section in the limit of D 2: Since the diagram behaves like (f)''^^, it becomes large when 
the dimension of the embedding space is large. When working in an e-expansion, d and D are 
related hy e = D — d{2 — D)/2 and sending £ ^ before D ^ 2 imposes that d ^ oo. Thus 
one would like to work at finite d, thus finite e. The diagrams are UV-convergent for d < dc- 
They are cut-off in the infrared by a finite membrane size. Alternatively, we can use a "soft" 
IR cut-off, that is we introduce a chemical potential r into the Hamiltonian and sum over all 
manifold sizes. For polymers (D = 1) this scheme corresponds to a "massive scheme" in an 
0(A^)-symmetric, scalar ^''-theory as can be seen via a de Gennes Laplace transformation on 
the = component limit. (The correlator in momentum-space reads p^) In the following 
we keep a hard IR regulator and calculate the two-loop diagram in fixed dimension. We use 
the analytical expression given in (4.16). For the numerical integration in 1 < D < 2 we 
proceed as before. The result for different embedding dimensions d is shown graphically on 
figure 5.1, where also some explicit values in the limit of D = 2 are given. Note, that the exact 
factorization-property in the limit of D = 1 is only valid in the MS-scheme or in a massive 
scheme with suitable soft cut-off as mentioned above. However, even the method is not tailored 
to capture that feature, the obtained value cu = 0.6 in the case of (D, d) = (1, 1) is not too far 
from the exact value 1/2. We expect that higher order terms would again reproduce the exact 
value of = 1/2. Also note that results for membranes (D = 2) are close to those from the 
e-expansion, especially for small d. 




d 




0.5 


-0.437 


1 


-0.410 


1.5 


-0.386 


2 


-0.358 


2.5 


-0.331 


3 


-0.299 


3.5 


-0.265 


4 


-0.229 



Figure 5.1: 2-loop diagram in fixed dimension for embedding dimensions d—Q.5 to d—A (from bottom to top). Tlie 
curves are only shown for values of D such that e>0. 



5.2 Self-avoiding membranes 



Up to now we were considering a D-dimensional polymerized "phantom" membrane interact- 
ing with a point-like impurity in rf-dimensional embedding space modeled through a (5-potential. 
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The infinite fractal dimension of the free manifold was causing problems, and since all opera- 
tors of the type (— Ar)"5(r), n>0, attain the same relevance in D = 2, it is not clear whether 
the model will break down. Physically, the problem is much better defined for self-avoiding 
membranes. From field-theoretical calculations, which have recently been refined up to two- 
loop [26,27], we know that two-dimensional self-avoiding manifolds embedded in three dimen- 
sional space have an anomalous dimension of u* ^ 0.85. We now try to overcome the problem 
of an infinite Hausdorff-dimension by a crude simplification: We approximate the self-avoiding 
manifold by a Gaussian theory with the same scaling dimension, that is a Hamiltonian of the 
type [16] 

= T^ [ lrix)i-A,)'/\{x) , (5.1) 



where 

r : X e — > r{x) E R'^ 
as before providing the two-point correlator in our normalizations: 

Cix,y) = \x-y\'~^ . (5.2) 

The scaling dimension is 

k-D 

V = . (5.3) 

We recover (3.9) by setting k = 2, but the model can be continued analytically to any real value 
of k, with k>2. Setting 

k* = 2u* + D (5.4) 

we get a Gaussian manifold with identical scaling dimension as obtained for self-avoiding crum- 
pled membranes at two-loop. The critical embedding dimension of the interaction then reads 

^ _ '2D _ D ^^^^ 

k*—D u* ' 

that is the interacting is relevant in d<dc. For membranes dc=-f^ ~ 2.4. All operators of the 
type (— A"5'^(r)) are at least naively irrelevant for n > 0, as 

[i-A^6\r))] = -u*{d + 2n) (5.6) 
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and the corresponding coupling has in inverse length units the dimension 

[^J")] ^D- v*{d + 2n) < , n > . (5.7) 

We may repeat the calculation in the massive scheme setting the free scaling dimension i/*=0.85 
and {D, d)—{2, 1). The result is shown in table (5.2). 

Let us stress that this method can not be turned into a systematic expansion, since perturbation 
theory neglects all effects of the non-Gaussian nature of the nontrivial crumpled state. We only 
know that for li — cxd the Gaussian variational approximation for the Hamiltonian, i.e. (5.1) 
becomes exact [33-35]. 





I 




d{9*) 


1-loop 


1.000 


1.15 


2.71 


2-loop 


-0.379 


1.49 


3.10 



Table 5.1: Two-loop results for self-avoiding membranes {D=2) approximated through a Gaussian theory having 
identical scaling dimension i/*=0.85. Results are shown for the corresponding diagram, the correction to scaUng 
exponent w and the contact exponent 6. 



6 Summation of the perturbation series in tlie limit D = 2 

In the last section we introduced a "massive" scheme allowing to perform the renormalization 
procedure in fixed dimension. As can be seen in figures 5.1 and 5.2, the limit D — > 2 can be 
taken in the 2-loop diagram and is smooth. Surprisingly, this limit can be calculated analyt- 
ically: Setting z/ = in (4.16), the integrand becomes constant, and thus trivial to calculate. 
This is a striking property. Even more, if we slightly modify the regularization prescription, 
this property holds to all loop orders in perturbation theory. Below, we will give an analytic 
expression for the complete perturbation series, which allows for analyzing the strong coupling 
limit. Furthermore, we show how corrections in (2 — D) can be incorporated. We will see that 
only the latter depend on the explicit cut-off procedure, a point which will be further clarified 
below. As a non-trivial check, when extrapolating back to polymers (D = 1), the corresponding 
result for uj is approximately reproduced. 



6.1 A^-loop order 

In order to calculate the iV-loop contribution to the perturbation theory we need the SDE of 
N + 1 contact interactions. The normal ordered product of the corresponding vertex operators 
reads 

Af+l 

gifeir(a;i)gife2r(a;2) . _ ^ikN+irixN+i) _ .g«(E^=i^ '=n''(^n)) . J~J g|fci*;jC'(a;i-Xj) (6 1) 



We choose a:;v+i as root in the renormalization procedure (for the precise definition, see below) 
and therefore /cat+i to integrate over in order to obtain the 5-distribution (or its derivative), onto 
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which to project. We thus shift 

N 

kN+1 — > kN+1 - ^ kj 
and rewrite the quadratic form in (6.1) as 

N+l N ^ N 



— kikjC{xi—Xj) — kM+ikjC{xN+\~Xj) + kikjC{xi—Xj) 

i,j=l j=l i,j=l 

N N ^ N 

y ] kN+ikjC{xN^i—Xj) — y ^ kikjC{xN+i—Xj) + — ^ ] kikjC{xi—Xj) 

j=l i,j=l i,j=l 



N N 

^ ^ U{XN+1-Xi) + (J{1 

2 

j=i i,j=i 

(6.2) 

(6.1) then becomes (up to subdominant terms^ involving spatial derivatives of r) 

:exp [ikN+irixN+i)]- 

N N 



X exp 



(6.3) 



and the integration over fcjv+i produces the 5-interaction and its derivatives; the latter come 
from the first term in (6.2). They are subdominant and thus will be dropped. We abbreviate the 
quadratic form in (6.2) as {Dij), where the matrix elements are 

= ^ [C{XN+1-Xi) + C{XN+1-Xj) - C{Xi-Xj)\ , iy^ j , 

Da = C{xN+i-Xi) . (6.4) 

Integrating out the momenta /ci, ... /cat yields the operator product expansion of + 1 contact 
interactions as 

6''{r{^,))6''{r{x2)) . . . 6''{r{x^+,)) = [det(A,)]-'/' ^'(r(a;^+i)) 

+ subdominant operators . (6.5) 

6.2 The limit D ^2 and (2-i:))-expansion 

Whereas (6.5) is still completely general, we now want to study the limit D ^ 2 as well as 
perturbations above it. We start by remarking that the matrix ( Aj ) may be expanded in powers 
of 2-D: 

- {2-D) 



n=0 



'Here we argue in terms of an OPE. Note that in order to make these arguments rigorous, one should ask 
of whether these terms affect the effective action at constant background-field. This is not the case, since then 
Vr = 0. 
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where 



D'^f = (In'' \{xN+i-Xi)iJ,\ + ln'^\{xN+i-Xj)n\ - In" \{xi-Xj)iJ,\) , i^j 

L'(r)=/.-^nn"|(x^+i-x,V|, 
11= 1/L. Then, 

det{D) = exp[Trln(Aj)] 



(6.7) 



= det(D(°)) ■ exp[Trln[l + ^ 



(2-D) 



n=l 



:p(0)]-l2^(n)]] ^ 



(6.8) 



where [D^°^] ^ denotes the inverse matrix to D^°\ Let us expand the OPE-coefficient in (6.5) 
up to second order in 2 — D. Expanding first the logarithm. 



n=l 



nl 



= (2-L')[L'(o)]-^L'W + i^-^ [[^(o)]-!^^^) _ ([^(o)]-i^(i))2j +c>((2-L')3) , 
and inserting into (6.8) one arrives at 

[det(L')]-''/2 = [det(L'(°))]-'^/2 [l - ^^"^^"^ Tr[L>W]-^L>W 

"^""-^"^^-^"^^^^ +0((2-D)3) 



[2-DYd 



Tr[L>(°)]-^L>(2) _ Tr([D(o)]-iD(^))2 - -Tr2[D(°)]-iD(^) 



(6.9) 



The zeroth order is [det (-0*^°^ )] In order to proceed we need the latter and the inverse matrix 
[£)(o)]-i. Let us first calculate the determinant. From (6.6) we have 



/I 



1 1 



-2i/ 



\ 1 



1 

TV 



1 1 

2 • • ■ 



\ - - 1 / 

\ 2 • • • 2 / 



V 



/1...0\ 



VI. ..ly 



-2i/ 



-(l + TVP), 



(6.10) 



1 . . . 1 



being the projector onto (1,1,...,1). Denoting by dim Im(P) the dimension of the space onto 
which P projects one has the general formula: 

det(l + aP) = (l + a)di^I^(P). 
Since here dim Im(P) = 1, we find 

-2.7V 1 + ^ 



2N 



(6.11) 



H. Pinnow and K.J. Wiese, Interacting Crumpled Manifolds 



27 



The inverse matrix of (1 + A^^P) can be inferred from the ansatz 

(1 + Arp)(i + 6P) = 1 . 

Since P is a projector P^ = P and this implies 

(N + b + Nb)F ^0 . 

Therefore, b = -N/{N + 1) and 

p(°)]-^ = 2(l-^p)/.^^ (6.12) 

To obtain the first order contribution in {2—D) from eq. (6.8) we only need [D^^^]~^ D^^\ which 
reads 

([D<»>l-DW)y = Uoff - ^ f A'!') . (6.13) 
The trace of (6.13) can easily be performed, with the result 

^r(lD^%^D<'>)^(^f:Dlp--^f:f:(l-S.,)DlA^.'■'. (6.14) 

V i=l i=l k=l / 

In order to factorize the integrals as much as possible, we found it convenient to modify the 
regularization prescription (3.21). In each order of perturbation theory we have to integrate the 
expression (6.5) over internal distances. These integrals have to be regularized in the infrared 
through an appropriate IR cut-off. In the preceding sections we have cut off the integrals such 
that all relative distances between internal points were smaller than L = 1/ fi. In the following 
we replace this by the demand that all distances from the fixed, but arbitrarily chosen point 
xn+1, the root of the subtraction prescription, are bounded by L. Other distances are not re- 
stricted (see fig. 6.1). Alternatively, one could calculate on a closed manifold of finite size. Of 
course, then the correlator has to be modified. Work is in progress to implement this procedure, 
which has the advantage of being more systematic [29] . 

To simplify the calculations, we further introduce the following notation: 



fi^ii T ■ ■ ■ 1 ^ik) '■ — / ' fi-^h ) • • • ) 

J XI Jxn 



(6.15) 



with new normalization 



1'-= f^J d^'xeUxlKl), (6.16) 



which is chosen such that the overbar can be thought of as an averaging procedure, and espe- 
cially 

T = 1 . (6.17) 
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Figure 6.1: Regularization scheme for the N-loop diagrams 



Thanks to our regularization prescription the integral of (6.14) over internal points can be re- 
placed by IjL~^^ (for the integration measure L^^) times 

2N' 



Tr([D(o)]-iD(i)) 



In I {xn+1 - Xi 



1 + N 

/2N{N -1 

~ V 1 + 



2N 

~ 1 + N 

Introducing a diagrammatic notation 



In \xn+i -Xi\ + 



In \xn+i - Xi\ - -\n\Xi - Xj 
N{N-1) 



1 + N 



1x1 I Xij^ X 7' I 7 • 



:= In I^AT+r 



and 




• — In I X j X i I ^ 



(6.18) 



(6.19) 



(6.20) 



the A/^-loop integral reads up to first order in 2— I? 



-d/2 p 
X 



2^ ' \ 1+N 



+ 0{{2-Df 




+ 



N{N-1) 
1+N V>j 




(6.21) 



We now turn to an explicit calculation of the diagrams. To first order in 2— D we have two 
contributions. The first (6.19) is 




/ In \xN+i-Xi\ ^ D 

Jxi Jo 

The second diagram (6.20) is 

jy2 rl 



^ dxx^ , 1 
X D 



--, as D ^ 2 . (6.22) 




JO 



2 I ^ X ^ d. X j In I X ^ X j I ^ 

Or 
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where we have put vectors to stress that Xi — xj is not a scalar. It suffices to evaluate this in 
D = 2. We find ^ 

= -T ■ (6.23) 




Combining (6.21), (6.22) and (6.23) we finally obtain for the renormalized coupling g up to first 
order in^ 2 — D: 



X 



1 - 



d / 2N 



2 \ 1+N 




N{N-1) 




+ 



-90 



■e\N 



1 + {2-D) 



1+N \X 
d f N N{N-1) 



This can also be written as 

-502'^/^-^)^ 



2 \1+N 4(l+iV) 



+ ■ 



(6.24) 



92'/' = - E 



Af=l 



TV! A^<^/2 



-(2-i^)- £ (z^2'^/V^)^(iV-l)(iV-2) 



A?"=l 



]\Jd/2+l 



(6.25) 



Absorbing a factor of 2*^/^ both in g and gfo and defining the dimensionless coupling z := g^jji ^, 
we obtain the final result 



9 



-E 



-Z] 



d 



N=l 



-^)^(iV-l 



2 ^ iV! 7V'^/2+ 

N=l 



Ar=i 



jY! J\[d/ 2+1 



+ 0{2-Df . 

This formula is the starting point for further analysis in the subsequent sections. 



(6.26) 



6.3 Asymptotic behavior of the series 

In the following we are interested in the limit of large z (strong repulsion), which also is the 

scaling behavior of infinitely large membranes. We therefore need an analytical expression for 
sums like (6.26) in the limit of large z. An example is the leading order expression 

°° i-z^^ 

Af=l 

^The expression in parenthesis is the inverse renormalization Zg-factor which relates the bare coupUng qq to 
the renormalized couphng g according to gZgfi^ = qq. 
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For the following analysis we need an integral representation of the above series. We claim that 
for all k,d> 



E 

Ar=o 



\N 



dr r^^^ ^ ex-p[-z e - kr] . 



(6.28) 



This can be proven as follows 

-j^ /"OO 



r(i) Jo 



^ N=0 

OO 

E 



/ drr'/^-'eM-ze-^-kr] = —rY.-l^ 



^d/2-l^-{N+k)r 



This integral-representation is not the most practical for our purpose. It is better to set r ^ s := 
e"'' which yields 



E 



,N 



N\{k + NY/^ r 



(2) -'o 



(6.29) 



This formula is already very useful for some purposes. It is still advantageous to make a second 
variable-transformation s ^ y := sz, yielding 



E 



(-^) 



AT 



"^mik + Ny/^ r(f)z^ 



In 



d/2-1 



Inz J 



(6.30) 



Finally we remark that we usually have the following combination 



N 



(Inz) 



d/2-1 



m{k + Ny/'^ 



dyy 



k-i 



1 - 



Iny 
Inz 



d/2-1 



(6.31) 



It satisfies the following simple recursion relation, which is helpful to calculate the /5-function: 

d 



(6.32) 



From (6.28) /fc( 2;) > for all /c, 0? > and the behavior for large z is obtained by expanding 
(1 -llf)'^'"' for small ^ 



(JnzY^ 

r(i) 



Jo 



dyy" 'e 



d 

- - 1 



Inz V2 



dyy'' ^ In y e 



+ 



{Inzy 



+ 0{e-') 



(6.33) 
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The result is 



_ (In zy/^-^m ( 



1 d-2V'{k) 



+ ... . 



(6.34) 



2, V Inz 2 V{k) 
For later reference, we note a list of useful formulas (7 = 0.577216 ... is the Euler-constant) 



(In zY'^-^ 7(ln zf/''-'' (67^ + tt^) (In zY'^ 



+ 



r(^) 



+ 



i2r(^) 



+ ... 



(In zY'^-^ (7 - 1) (In ^)'^/2-2 (g^(^ _ 2) + 7^2) (In zY''"''' 



+ 



r(^) 



+ 



+ 



2(in^)d/2-i (27-3)(ln^)'='/2-2 (67^ - I87 + 6 + 7r2)(ln^)'='/2-3 



+ 



+ 



r(|) ■ r(^) ■ 6r(^) 



(6.35) 

(6.36) 

+ ... . 
(6.37) 



Note that with the above notations, the relation (6.26) expressing gr as a function of z becomes 

g{z) = ft^\z) - (2 - D)pt\z) + (2 - D)pt\^) + 0(2 - Df . (6.38) 

In the next two sections we shall present two formalisms which give results for uj at next 
to leading order. It will turn out that the results catch qualitatively correctly the cross-over to 
D = 1; however they differ in the numerical values. This is a reflection of the fact that the 
(2— -expansion is not systematic. We expect that working on a manifold of finite size and 
using the exact correlation function on this manifold would lead to a systematic expansion. 
Work is in progress [29] to check this hypothesis. 

6.4 The RG-functions in the bare coupling 

The renormalization group /3-function is defined as 



9 



(6.39) 



Since it is hard to invert the function g{z), it is advantageous to study /3 as function of the bare 
coupling z 

Pig) = = -ez—g . 



To first order in 2 — £), we find from (6.38) using (6.32) 

m = -e 



ftiz) + {2-D){ + ) + 0(2 - Df 



Inserting the asymptotic expansions from the last section, this becomes 



(6.40) 



(6.41) 



r(i) 



{lnzf'-'-'^—^{lnzf' + ... 



(6.42) 
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-60 20 40 60 80 100 -2o 20 40 60 80 100 



z z 

Figure 6.2: /3-function (left) and w (right) as functions of the dimensionless bare couphng z for different dimen- 
sions d: d = 0.25, 0.5, 0.75, . . . , 3 (from top to bottom, respectively). Note, that /3 always has a fixed point at 
z=0 with uj{0)= — e. Furthermore, there is a fixed point for 0<d<2 at z=oo with a;(oo)=0, the latter remaining 
true for d>2. In d=2 (3 tends to —e as z ^ 00. Above, /3 diverges. 



where omitted terms have either one more power of (2 — £>) or 1/ In z. As shown in section 2, 
universal properties are given by the correction to scaling exponent uj, reading 

Let us first study D — 2. P{z) and uj{z) are plotted on figure 6.2. We see that for d < 2 the 

/3-function becomes for z ^ 00. For d > 2, P{z) has no zero, but the flow of g is still to 
infinity. In both cases u is given by c<j = Iim2;^ooa;(^). Taking only the leading term in the 
(2 — £))-expansion of /3{z), lj{z) is given by 

z^f^(z) d-2 

uj(z) = -E '^l,] \ ' = — > . (6.44) 

^ ' ft{z) 2hiz 

Let us now take into account the first order in 2 — D. Then the fixed point of (6.42) is at the 
finite value 

In.- = ^ . (6.45) 

Next, when trying to calculate uj{z), we face the following problem: Since we truncated the 
series for g at order one in 2 — L>, I3'{z) does not vanish at the fixed point, i.e. the zero of I3{z). 
This might lead to the conclusion that oj{z*) is always 00, an absurd result. However since 
this is a consequence of the truncation of the series, we follow the strategy also to expand the 
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denominator ofu{z) in powers of 2 — D. From (6.41) and (6.43) we obtain: 



u;{z) — —e 



+ {2-D\ 



+ {2-D) 
Inserting the asymptotic series, we find 



(6.46) 



d-2 2-D 

CUiz) — —£— h £ 

^ ' 2\nz 2 



(6.47) 



Inserting the fixed point (6.45), we arrive at 



u{z*)^e{2-D) 



4-d 



This should be checked against the exact result in D = 1, which reads 



(6.48) 



(6.49) 



Thus for d — > our resummation-procedure gives the exact and for d > an approximative 
result. Before analyzing the validity of the procedure, let us turn to the second scheme, namely 
the calculation in the renormalized coupling. 

6.5 Calculation in the renormalized coupling 

We start from g(z) given in (6.38). For large z this can be approximated by 



r(|+i) 



d {\nzr/'^' 

^ U r(f+2) 



(6.50) 



where we only retained the leading term from each series. Also note that this formula is only 
valid for d < 2. For d > 2 additional additive constants have to be added. Using (6.50), we can 
write In 2; as a function of g. To first order m2 — D, this reads 



lnz = ^2/<i^^_^~4/<i 



2 + d^ ' + 

We now write the /^-function in terms of g. Starting from (6.42) we obtain 

£ 



(6.51) 



P{~9) 



^ 2+d ^ 



(6.52) 



Fixed points are at the zeros of the /3-function. The non-trivial one is at 



9 



2+d 
_2(2-D)_ 



d/2 



(6.53) 
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Figure 6.3: /3— function in terms of the renormalized coupling g truncated at order 160, Pade-resummed, and 
plotting only that part for which the truncated series converges. (This can e.g. be tested by taking away the last few 
terms of the series.) This is compared to the asymptotic behavior (6.52) (proportional Xol/g for large g). d is set 
to 1, and we used the diagonal (80,80)-Pade approximant, which was find to converge best. (The non-resummed 
expression starts to diverge already at 5 « 1.8 at this order.) 



The correction to scaling exponent uj is simply obtained by evaluating the derivative of (3 with 
respect to g at the fixed point. In terms of g this is 



dg dg dg ^ dg 

d ^\ f,_d{2-D) 



2 J " 2+d 
Inserting the fixed point g* from (6.53) we find 



(6.54) 



uj{~g*)=e'^^^J^ . (6.55) 

Again, this is quite close to the exact result uj{D=l) = e in D = 1, and in effect exact in D = 1 
and d = 0. 

Let us finally point out that in the limit D = 2 the true asymptotic behavior of the /3— function 
in terms of the renormalized coupling g is obtained from the completely summed series (6.27) 
leading to (6.52) for large g. Conversely, if one tries to invert (6.27) and truncates it taking only 
a finite number of orders into account, it is at least possible to reach the asymptotic regime - 
however, for large enough g the truncated /3— function wildly oscillates and thus strongly de- 
viates from the true behavior. In figure (6.3) the Pade-resummed truncated /3— function up to 
order (7^^° in (i = 1 is compared with the exact, asymptotic flow-function. One notices that the 
truncated /9-function eventhough improved throuh a Pade-Resummation hardly gets into touch 
with the asymptotic regime. The same applies to the slope-function uj{g), which is not shown 
in fig. (6.3). 

Note that the above arguments suggest quite intuitively the behavior of the exact /3-function 
in 1 < < 2. Whereas iox D = 1, the /^-function is a parabola, and for D = 2 it decays like 



H. Pinnow and K.J. Wiese, Interacting Crumpled Manifolds 



35 




Figure 6.4: Qualitative behavior for the /3-function in D = 1, D = 2 and result anticipated for D fa 1.5. 

a power-law for large g at least as long as d < 2, the /3-function for values of D between these 
two extremes should cut the axes /3{g) = at a finite value of g, which for D — > 2 wanders off 
to infinity, thus by continuity forcing the exponent a; to go to for Z) ^ 2. 

6.6 An instructive case: d = 

Setting d = in the leading order series (6.27) allows for a simple analytic expression for g as 
a function of z and its inversion 



g{z) = 1 - e ^ 
z{g)^-Hi-g) 



(6.56) 
(6.57) 



This gives 



P{g) = £(1 - g) ln(l - g) 
a;(^) = -£(l + ln(l-^)). 



(6.58) 
(6.59) 



Fixed points are g* = and g* = 1 with (jj{g*=Q) = —e and u){g*=l) = oo. Note that the 
latter result is an artifact of d = 0. 



7 Discussion and conclusion 

In this article, we have discussed the perturbation expansion of a D-dimensional elastic man- 
ifold interacting via a ^-interaction with a fixed point in embedding space. This calculation 
can be done exactly for D = 1, but becomes non-trivial for D ^ 1. Interestingly and quite 
surprisingly, it simplifies considerably in the limit of D 2, if one decides to work at finite 
e. In that limit we were able to obtain an explicit expression for the renormalized coupling 
as function of the bare one, in terms of a non-trivial series. Analysis of this series shows that 
the fixed point lies at infinity in the bare coupling, a limit in which we were able to derive an 
asymptotic expansion for the renormalized coupling as a function of the bare one. This yields 
a vanishing exponent tu in the limit of D — > 2. Also it is important to note that this result is 
completely independent of the regularization procedure. This does no longer hold true beyond 
the leading order, which should be accessible to an expansion in {2 — D). We here constructed 
its first order in a specific regularization scheme. While this reproduces qualitatively correctly 



36 



H.A. Pinnow and K.J. Wiese, Interacting Crumpled Manifolds 



the known results in D = 1 (and even exactly for d — 0), it also shows by its dependence on the 
renormalization scheme (here working with the /3-function either in terms of the renormalized 
or bare coupling) that this expansion in (2 — D) is not systematic. Even though experience with 
this new kind of expansion is still lacking, this is likely to be caused by the use of a hard cutoff 
in position space while working with the correlator of an infinite membrane. It seems that only 
in an £-expansion this procedure is systematic. Work is in progress [29] to study the model on a 
sphere or torus of finite size, such that no further infrared cutoff is necessary. This should lead 
to exact results beyond the leading order and provide for the crossover between polymers and 
membranes. 

While results for the pinning problem are interesting in its own, the main motivation is 
certainly to obtain a better understanding of self-avoiding polymerized membranes. This model 
is well-behaved physically, since its fractal dimension is bounded by the dimension of the 
embedding space df < d. Preliminary studies [29] indicate that this problem can also be 
attacked by the methods developed in this article. This would be very welcome to check the 
2-loop calculations [27,26] and the large-order behavior [36] on one side and numerical results 
(e.g. [28]) on the other. 



A Universal l/r - repulsion law for polymers 

We have shown in section 2.2 by a scaling argument that as long as r ^ L", the restricted 
partition function of a manifold pinned at one of its internal points scales as Zoo{r/L^) ~ 
[r/L^y, with the contact-exponent given by (2.25). It is instructive to prove this for the special 
case of polymers m d = 1. According to (2.1) and (2.6) this corresponds to a polymer in 
3d-space interacting with a (5-potential on a plane. 

Let us start from the (bulk) polymer propagator being defined as the conditional probability 
of finding the internal point s' at the position r/ when starting in s at r^. The propagator can be 
written as a functional integral for a restricted partition function of the free chain according to 



ZsAnVs) = j V[r] S{ris)-ri)d{r{s')-rf) e'^M 

Jk 



'k 

The propagator (A.l) posesses the Markov property 



Zs,s"{ri\rf) ^ J d'^Vm 2s,s'{ri\rm) 2s',s"{rm\rf) . (A.2) 

In the following we are interested in the restricted partition function (A.l) in presence of a 
short-range interaction modeled through a 5-potential, which is situated at the origin of the 
embedding space. The restricted partition function in the presence of the interaction will be 
denoted by g,(rj|r/). We consider the case, where the chain is fixed at its ends, such that s=0 
and s'—L. We furthermore switch to a grand-canonical ensemble and denote 

POO 

Z?{n\rf):= dLZl\{n\rf)e--'^ , (A.3) 

^0 
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T being some chemical potential. Zf^{ri\rf) can be expanded in a perturbation series in the 
coupling qq according to 



n=0 



TheZf ■'*(rj|r/) are 

oo 

''^(r,|r/) := j dL j Zq,,, (r,|0) Z,,,,,(0|0) 



(A.4) 



X . . . 



0<xi<...<x„<L 



-tL 



- n-1 






11 



...xZ,„_„,„(0|0)Z,„,L(0|r;)e 



= Z°(r,|0)Z°(0|0r^ZO(0|r^), (A.5) 
where we have used that after integration over L, the integrals factorize; we further have defined 

poo 

Jo 



poo p p ^ikirt—ri) 

= / ds e-ik^+r)s+ik{rf-ri) ^ / ^ _ 

Jo Jk Jk k"^ + T 

In d—1 the integration over the momenta k in (A.6) leads to 



(A.6) 



(A.7) 



Let us now lock at what terms appear in the n-th order coefficient in (A.5): From the initial point 
of the chain a factor Zr{ri\0) is contributed, while the final point is coming up with Z^(0|r/). 
Furthermore, the internal points give n—l powers of Z^(0|0), where 



2r{0\0) 



The (n=0) -coefficient is special: Clearly, one has 

Zf'°(ri|r/)=Z,(ri|r/). 
Inserting (A.8)-(A.9) into (A.5) we obtain for (A.4): 



(A.8) 



(A.9) 



(A.10) 
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The limit of ^ can be taken and leaves us with the beautiful result 

Z^{ri\rf) = (e-V^lrf-nl _ Q~V^{\rf\+\ri\)^ . (A.ll) 

This is the propagator of a scalar field with Dirichlet boundary conditions. It has applications 
when studying critical phenomena of e.g. an Ising magnet in half-space [37,38]. Performing the 
inverse Laplace-transformation of (A.ll) we arrive at 

^rinlrf) = jjj^ ie-^w- - e j . (A.12) 

We conclude that a 6- potential acting in a hyper-plane suppresses all configurations which 
penetrate or touch it, when taking its amplitude to infinity. 

Finally, in order to prove the universal repulsion law, we start from (A. 12) and integrate over all 
final positions r/: 

oo 

2r{r^) = J drj Zr{n\rf) = 2^ (^1 - erf , (A.13) 

— oo 

erf denoting the error function. For small arguments we have 

l-erf(x) ~ X. (A. 14) 

Furthermore, Z^(rj) denotes the restricted partition function of a chain pinned at one of its 
ends in r^. To obtain Zoo{r / W) in (2.23) we have to evaluate: 

Z^ir/n -ll'ds Zr{r)Zr_Ar) - (^) (A.15) 

according to (A. 14) in the scaling regime r^L^/^. We thus find that 9 = 2, which confirms 
(2.25) for Z) = 1 and d = 1. The above proof can be extended to arbitrary dimensions Q<d<2, 
where s>Q for polymers. Then, ^ = 4 — 2d. 

It would be nice to make the same arguments for membranes. However, the proof is based 
on a drastic simplification, which only occurs in D=l and shows up in the factorizability of 
loop diagrams as in (A.5). This has no extension to manifolds of internal dimension D>1. 



B Conformal mapping of the sectors 

In order to calculate the two-loop diagram efficiently, one wants to write it as an integral over a 
finite domain only. To do so we need the technique of conformal mapping of the sectors, which 
also serves for analytically continuing the measure of integration to internal dimensions D < 1. 
This technique has been extensively used and well documented in the context of self-avoiding 
tethered membranes [30,27,4], but we repeat the presentation here for completeness. 



H. Pinnow and K.J. Wiese, Interacting Crumpled Manifolds 



39 




Figure B.l: The sectors A, B and C 

Generally, in evaluating the three-point divergences we need to integrate over some domain in 
the upper half -plane (see 4. 11). The measure of integration reads 



oo 



/ = ^ / I dy, {y,r^' (B.l) 



-oo 



and the integrand is a function / of the three distances a, h and c between the points O — (0, 0), 
H = (— L, 0) and Y = [yi, 7/2), given explicitly by (see fig. B.l) 

a = L = fixed, h = + {y2Y , c = ^/(L=yO^ + Ty^ ■ (B.2) 

In our problem / is homogeneous, of degree A, but not necessarily symmetric: 

/(«;a, K,b, Kc) — K~^f{a, h, c) . (B.3) 

Let us now explicitly show, how sectors can be mapped onto each other. Specializing to the 
mapping i3 — .4 we choose a coordinate system as given in figure (B.l). The mapping B — > 
A is mediated by a special conformal transformation, the inversion with respect to the circle 
around O and radius L. In complex coordinates this is (Y is the complex conjugate of Y) 

7-2 7-2 ^ 7-2 

Y^Y^^^Y- ^Y^Y^ (B.4) 



such that 



62 / X^^^/ \ p 

This change of coordinates gives a Jacobian for the measure (B.l) 



2N D 



(B.6) 
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and 



It can easily be seen that the mapping B ^ Ais one to one: First, be Y={yi, 1/2) G B, that is 
b > L and b > c. Then, 

and 

b' = L'^ < . (B.8) 

Since 6^, < L^, F e ^. Second, for the inverse mapping, Y{Y e A) G <B is checked as 
follows: 

^ b'' + L^- 2y,L ^ b'^ + L''- 2yS < 6' , (B.9) 

since 

r3 

= 6^ - 2yiL + < ^ - 2yi^ < 
7-2 

6^ = > . (B.IO) 

Since b > L and 6 > c, F e S. Finally, let us look at how the whole integral transforms. 
First inserting the transformed variables including the Jacobian (B.6) and second using the 
homogeneity (B.3), we arrive at 

l^jj fl^L,j,jj= J l^jj f{b,~a = L,d). (B.U) 

Let us restate the above calculation in a completely geometric interpretation as sketched in the 
figure. The mapping consists of two steps: 

• The rescaling with respect to O by a factor L/b which maps Y onto Y' and H onto H'. 
(B.3) implies that / is changed by a factor (-|) 

• A mirror operation, which maps Y' onto H and H' onto Y, leaving invariant the origin 
O. This operation is a permutation of the first two arguments of /. 

The mapping Y Y corresponds to the special conformal transformation (B.4). 
Analogously, we find that 

f{a = L,b,c)= j f{c,b,~a = L). (B.12) 

y^c YeA 
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